Rheology of Active Filament Solutions 
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We study the viscoelasticity of an active solution of polar biofilaments and motor proteins. Us- 
ing a molecular model, we derive the constitutive equations for the stress tensor in the isotropic 
phase and in phases with liquid crystalline order. The stress relaxation in the various phases is dis- 
cussed. Contractile activity is responsible for a spectacular difference in the viscoelastic properties 
on opposite sides of the order-disorder transition. 
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Soft active systems are exciting examples of a new 
type of condensed matter where stored energy is contin- 
uously transformed into mechanical work at microscopic 
length scales. A realization of this are polar filaments 
interacting with associated molecular motors in the cell 
cytoskeletoniS. These systems are characterized by a 
variety of dynamic and stationary states which the cell 
accesses as part of its cycl o 3 i 4 i 5 . Recent experimental and 
theoretical studies of the dynamics of solutions of active 
filaments have focused on the formation of both homoge- 
neous and inhomogeneous states with spatial structures, 
such as bundles, vortices or asters^&L&SiiSiiiiiSii&ii. 

In this letter we study the effect of motor activity on 
the rheological properties of active solutions under an 
externally imposed stress. Understanding the viscoelas- 
ticity of cells, and cellular extracts in-vitro, is a very 
important problem currently receiving a lot of experi- 
mental attentioni^*iSiiLi2iiS. From a theoretical point of 
view describing the mechanics of the cytoskclcton in its 
full complexity remains very challenging. As a first step 
in this direction we focus here on using methods from 
polymer physics to understand the effect of motor ac- 
tivity on the viscoelasticity of a dilute solution of long 
stiff biopolymers. A phenomenological description of the 
rheology of isotropic suspensions of active particles near 
the isotropic-nematic transition was proposed recently by 
Hatwalne et alw The present work provides a micro- 
scopic basis for their results and a general framework for 
analyzing the viscoelastic behavior of active solutions in 
both isotropic and ordered states. 

Our model makes several testable predictions. First we 
find that, as suggested in RefiSS, activity yields a con- 
tribution to the viscosity of an isotropic solution that 
diverges at the isotropic-nematic transition (see Fig. QJ. 
Such a diverging viscosity is reminiscent of an equilibrium 
liquid-solid transition rather than a liquid-liquid transi- 
tion, and is a direct consequence of the active-stresses. 
Unlike a liquid-solid transition, however, the divergence 
is localized at the IN transition and the viscosity is finite 
in the nematic phase. A second novel signature of activ- 
ity in the stress relaxation is found in the nematic phase, 
where the shear stress acquires a nonequilibrium contri- 
bution proportional to ATP (Adenosine Tri-Phosphatc) 
consumption rate that remains finite for zero deformation 
rate. In other words, an active nematic solution is driven 
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FIG. 1: (a) The steady-state shear viscosity on both sides 
of the IN transition for passive and active nematics.(b) The 
'motility' parameter, \ ( see Eq. 1201 ') relating the magnitude 
of the active component of shear stress to the 'activity' n ( 
~ ATP consumption) Inset: phase diagram showing the IN 
transition. 



into a state with a non-vanishing macroscopic stress by 
the energy input from ATP hydrolysis, even in the ab- 
sence of an externally applied mechanical deformation. 
This non-equilibrium contribution is also present in the 
normal stresses in the nematic phase. 

We consider a suspension of polar filaments in a thin 
film of height comparable to the length of the filaments 
(quasi-two dimensions) and a constant density m of mo- 
tor clusters. A concentration p(r, t) of filaments is sus- 
pended in an incompressible solvent of viscosity r/o char- 
acterized by a fluid velocity v(r, t), with V • v = 0. Mo- 
mentum conservation yields 
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with p s the (constant) density of the solution. The sol- 
vent contribution to the stress tensor is 

V • cr s = j) V 2 v - Vp , 

withp the pressure. The filament contribution, cr* , must 
be determined in terms of the driving forces (velocity 
gradients, Kij = djVi, and motor activity, p, ~ ATP con- 
sumption rate) and the conserved and order parameter 
fields describing the filaments, 



a s = f {p,v,Sif,K,p) , 
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where p and Sij arc the local polarization and nematic 
alignment tensor that describe the oricntational order of 
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the filaments. The derivation of this constitutive equa- 
tion from a model of filament dynamics is one of the 
central outcomes of work. 

The filaments are modeled as hard rods of fixed length 
I and diameter a (I 3> a) at position r with filament 
polarity characterized by a unit vector u. The filament 
contribution to the stress tensor is^i 



V 
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where T(s) is the hydrodynamic force per unit length 



r l/2 



ds. 



on a rod at position s along the rod, (..) s = J_^ 2 ' 
and c(r, u, t) is the concentration of polar filaments with 
position/orientation {r, u}. The force on the rod is spec- 
ified by its interaction with the solvent, other rods, and 
the motor clusters. For low Reynolds numbers Re <C 1, 
viscous effects dominate inertia and we can set the left 
hand side of Eq. Q to zero. 

We calculate the force per unit length ^(s) by de- 
composing a rod into a sequence of beads of diameter 
a and solving self-consistently for the flow field around 
the rodSI on scales much bigger than a. The stress due 
to the filaments is (to C(V 2 )) 



V-*f(T,t) = /f(r,u,t) 
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where 
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with U x (r, u) = k B T J_, J c(r + £, u') the excluded vol- 
ume potential. The force density has contributions from 
fluctuations or diffusion (both thermal and non-thermal 
- hence the active temperature T a ^ T ), excluded vol- 
ume, and motor activity. There is also a viscous contri- 
bution to the stress proportional to the velocity gradient 

(Cx=47rtj ? l/ln(i/o))- 

The active force and torque are given by 



F a (r,u 1 ) = -m / (C(ui)-v (l;2)c(r + £,u 2 )) , 

J u 2 

T„(r, fix) = -m f <CrW a (l; 2)c(r + £, fi 2 )} sliS2 , (5) 

J U 2 



where mo 



\ (1;2) = (si,ui;s 2 ,u 2 ), £ = uisi - 



U2S2 and £(u) = C-L — flu) + C|| uu ) with Oi and 
£ r friction coefficients. The angular velocity is taken 
as u) a = 2 [70 + (fii • 62)71] (fli x "2), with 70 and 71 
motor-induced rotation rates proportional to ATP con- 
sumption, and tends to align filaments^ The motor- 
induced translational velocity has been derived from a 



model of motors walking along the filaments at a mean 



rate /?. It has the form v a (l; 2) 



V m , with 
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(u 2 - fii) 
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where a = a(l + fii • u 2 ) and = (3(1 + Ui • fi 2 ). The 
parameter a ~ f3(a/l) is controlled by spatial inhomo- 
geneities in the motor stepping rate. Momentum con- 
servation yields expressions for A and B. For long thin 
rods with (j_ = 2£|| = 2f, to leading order in Ui • u 2 , we 
find A = -\p-a(s 1 + s 2 )/2]/12 and B = a(si - s 2 )/24. 
When evaluating the contribution to the stress tensor, 
only terms up to first order in Ui • u 2 are retained in the 
active force C( u i) • v a(l; 2) exerted by a motor cluster on 
the filament in the first of Eqs. Ip|l. This approximation 
only affects the numerical values of the coefficients in the 
stress tensor, not its general form. 

The concentration c(r, u, t) of polar filaments satisfies 
a local conservation law, 

dtc+V -(vc)+TZ-(ljc) + V ■ 3 + TZ- J = 0, (6) 

where uj = u x k ■ u and 1Z = u x Jt. The transla- 
tional and rotational currents in Eq. © contain diffu- 
sive, excluded volume and active contributions^ J = 
— C _1 (u) -f(r, u) , J = — C J T 1 T"(r, fi) , with f and r given 
by Eqs. J3J. 

The conserved and broken symmetry fields 
arc the density p(r,t), polarization p(r, t), and 
nematic order SV/(r, t), defined as moments of the proba- 
bility distribution, J. c(r, u, t) = p(r, t) , J. uc(r, u, t) = 

p(r,t)p(r,t) , /. Qijc(r,u,t) = p{r,t)S lj {r,t) , where 
Qij = Uiiij — Continuum equations for these 

fields are obtained from Eq. © by the coarse-graining 
procedure described in RefiSi. 

The constitutive equation for the stress tensor cr? is 
obtained by evaluating the right hand side of Eq. For 
simplicity we consider spatially homogeneous solutions in 
the presence of a constant velocity gradient . To lowest 
order in gradients, the three contributions to the stress 
tensor of the filaments are 



4 (r, *) = 0% (r, t) + at (r, t) + (r, t) , 



(7) 



with passive and active contributions, <x p and cr A , given 
by 

cr[j = 2k B T a p (l — S )Sij--^—(piPj-\p 2 5 i j) 

l\ Pin' APtp I 
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pk B T a p 2 
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where p 



i8D 



Pip = D r /(m 7oZ 2 ), and p IN 



P7v/[l + P7v' 2 «^o7i/(4£>r)] is the density for the isotropic- 
nematic (IN) transition at finite motor density^, with 
Pn = 37r/(2/ 2 ) the density of the IN transition in passive 
systems. Finally, the viscous contribution to the stress is 
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where tq = l/D r and D r = ksT a /( r the rotational diffu- 
sion constant. The time scale ta diverges as we approach 
the active IN transition. In a sheared sample, the total 
stress (filaments + solution) is 

°xy = fje+2k B T (1 — ) pS xy +^-pk B T p 2 S xy , (16) 

\ Pin J » 

(10) where fj = ?/o(l H — )• The first contribution is from 



with K?j = (ttij + Kji)/2. 

In agreement with Ref^i, we find that active units 
generate force dipoles in the fluid yielding contributions 
to the stress which are equilibrium- like, i.e., have the 
same form as those appearing in an equilibrium solution, 
but with new contractile stresses (a > 0) which have no 
analogue in their equilibrium counterparts. 

For a homogeneous solution p = constant and the 
equations for polarization and nematic order parameter 
are obtained by averaging Eq. © over u, 



d t pi 



dtSij 



Pi 
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with pCli = J duuiTt-{bjf) and pClij = J duQijTZ- (cj/) 
We find 



-KijPj + - 



K i]Pj + 2 K kkPi 



(13) 



Clij — — ^ij K kk\ { K ikSkj + KjkSki*) 



+ - 



SijKkk + SijSkinu + 2(SikK S kj + Sjkn s ki ) .(14) 



The homogeneous (bulk) steady-states are obtained 
by setting the rhs of the order parameter equations to 
zero. Wc find (I)sotropic (p = 0, Sij = 0), (N)ematic 
(p = 0, Sij ^ 0) and (P)olarizcd (p ^ 0, Sij ^ 0) phases. 
It can be shown using Eqs. 1)819(1 that in a stationary bulk 
fluid the passive contribution to the stress tensor is iden- 
tically zero in each phase, while the active contribution 
is non-zero. In the following we consider the geometry 
of pure shear flow with = i8i X 5j V and discuss the 
linear viscoelastic response of the active solution in the 
isotropic (I), polarized (P) and nematic (N) phases. 

a. Isotropic phase In the isotropic phase close to the 
IN transition a shear flow builds up nematic order, yield- 
ing a non-zero value for Sij to 0(e). The shear stress is 
linear in the strain rate, a xy = 0(e), while the normal 
stress is quadratic, a xx — a yy = O (e 2 ) . To linear order 
we obtain (dt + 1/ta)S xv = e/4, with 
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the solvent and the viscous stress, the second is from the 
passive stress, and the third from active stresses. 

For an oscillatory applied shear we define the frequency 
dependent shear viscosity o~ xy (u>) = r)A(to)e, with low fre- 
quency limit 



Va 



~^k B T P( 16 



(17) 



This diverges as we approach the I — N transition (see 
Fig. 2]). This behavior should be contrasted to that 
of a passive solution of rods near the equilibrium I- 
N transition. In this case the stress relaxation time 
Tp = D^ 1 /(I — p/ pn) also diverges at the IN transition, 
but the zero- frequency viscosity, rj = fi+ksT p/ (&D r ), re- 
mains finite, as required by the Fluctuation-Dissipation 
Theorem. 

b. Nematic phase In the nematic phase, there is 
the possibility of alignment of the nematic director by 
the shear flow. The director n and the magnitude S 
of the nematic order parameter are defined by Sij = 
S (hihj — /2). The equation of motion for the director 
in a steady shear flow is obtained from Eq. Q12|). 



dtfii = {Sij - h t hj) n k 



s 



(18) 



Defining n = (cos 9, sin 9), we obtain a steady state so- 
lution for the director orientation given by cos 29 = 2S. 
For S > 1/2 there are no steady-state solutions, possibly 
pointing to the existence of periodic or chaotic solutions 
characterized by "wagging" or "kayaking" of the nematic 
director—. 

The steady-state expression for the stress tensor in the 
active nematic state is obtained from Eqs. JSJl and as 



k B Tp { Vkj_ 
2D r V P 



(19) 



from which the 6 Leslie coefficients of the active nematic 
can be obtained using Eqs. lf§)) . (fTT)|) and frf f 21 ^ . 

The novel nature of the constitutive equation is best 
illustrated by the simple case of an active nematic in the 
flow-aligning regime (S < 1/2), where the steady-state 
shear stress is given by 



o xy = m + xn , 



(20) 



where r, = fj + ^Lp (± - S 2 ) 1 + D r r' A ^\ and X = 



lk B Tp 2 SJ\ - S 2 , with 1/t' a = 8D r (p/p IN - 1). The 
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magnitude of nematic order S relaxes on a time scale 
t' a while the director relaxes on the shear time scale e 
leading to non-monotonic stress relaxation. 

The signature of this constitutive equation is a shear 
stress that does not vanish for zero deformation rate. 
This is because the active filament system is being driven 
out of equilibrium by two sources of energy. One exter- 
nal, due to the shear, and the other internal, due to the 
activity of the motors (see Fig-QJ. The viscosity diverges 
as the IN transition is approached from the nematic side, 
but it decreases dramatically as one goes deeper into the 
nematic phase (see Fig.^). This is a direct consequence 
of flow alignment as S increases with density. All cor- 
rections to r\ vanish at S = 1/2, where the flow-aligned 
regime ceases to exist. 

Another signature of a nematic phase is that the nor- 
mal stresses are of first order in the shear rate. In ad- 
dition, in an active nematic one obtains an anomalous 
constitutive equation, 

Oyy - <T XX = rj N e + XNfJ- , (21) 

where m = %£ P SyJ± - S~ 2 and X n = lk B Tp 2 2S 2 . 

Both the anomalous stresses (shear and normal) should 
be easily observed by performing linear rhcological ex- 
periments on active filaments in the nematic phase at 
varying shear rates. 

c. Polarized phase In a polarized state, with p = 
PoPi a uniform density and a constant velocity gradient, 
the unit vector p satisfies the dynamical equation 

dtPi + UijPj = XSfjKjkpk , (22) 



where A = 1/2, = (n Jl -K ij )/2 and Sfj = (6ij — piPj). 
Voituriez et al. — have recently suggested that sponta- 
neous flow can be obtained in active polar materials, cor- 
responding to a state with finite velocity gradients and 
flow alignment. Flow alignment, however, can only oc- 
cur for A > 1, while the present (low density) calculation 
yields A = 1/2. Corrections to A can be obtained by in- 
corporating the fact that fact that the alignment tensor 
is slaved to the polarization field. However an analysis 
of such corrections shows that they fail to increase A to 
values larger than 1, suggesting that no steady uniformly 
flowing polarized state exists for a thin film in the low 
density limit presented here. An interesting alternative 
is the possibility of periodic or chaotic states^. 



In summary, we have used a molecular model to study 
the macroscopic mechanical response of active filament 
solutions in both isotropic and ordered states. Motor 
activity leads to a novel coupling of mechanical proper- 
ties to order and to anomalous constitutive equations in 
the liquid crystalline state. The theoretical framework 
developed here can be generalized to consider stress in- 
homogeneities. This is relevant for understanding the 
microrheology of active filament systems where new be- 
havior is expected even in the isotropic rcgims2&. 
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